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Rastall generalized Einstein's field equations relaxing the Einstein's assump- 
tion that the covariant divergence of the energy-momentum tensor should vanish. 
His field equations contain a free parameter a and in an empty space, i.e. if = 0, 
they reduce to the Einstein's equations of standard general relativity. We calculate 
the elements of the metric tensor given by Rastall' equations for different a assum- 
ing that T^ u to be that of a perfect fluid and analyse these model solutions from 
the point of view of Mach's principle. Since the source terms in Rastall's modified 
CN ■ gravity equations include the common energy-momentum tensor T^ u as well as the 

expressions of the form (1 — a) g^T/2, these source terms depend on metric deter- 
mined by the mass and momentum distribution of the external space. Likewise, in 
the classical limit, the source term of the corresponding Poisson equation for a / 1 
depends on the gravitational potential in the sense of Mach's principle. 
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I. INTRODUCTION 



q-i General relativity assumes that particle rest masses are space-time constants independent 

jLi I of the space-time metric in which they occur. The experimental evidence relating to this 
assumption is fragmentary and far from establishing it on a firm basis. The space-time con- 
stancy of particle mass is a typical feature also of Newton's mechanics. Mach in his critique 
of the Newtonian mechanics [2| pointed out to the fact that mass and momentum distribu- 
tions of the external space should be linked with inertia of a mass body. This approach to 
treat problems of mechanics is commonly called as Mach's principle |2(. Many versions of 
Mach's principle can be found in the literature (see, e.g. 0,0,0]) and unique, satisfactory 
"canonical" formulation of the principle does not seem to exist as yet. Nevertheless, there 
seem to be a fair consensus about the consequences which should follow from it. Mach's 
principle requires that the presence and motions of large external masses in the universe 
must have some definite effect upon the mass objects occurring in it. As is well-known, 
Mach's principle is not incorporated in the Einstein equations of classical general relativity. 
To be included, the standard general relativity equations have to be modified. In the liter- 
ature, there are several attempts to modify the Einstein equations in such a way that they 
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incorporate at least some requirements following from Mach's principle (see, e.g. |lOj). 

One of the fundamental assumptions of the Einstein theory of gravitation is vanishing 
divergence of the energy-momentum tensor 

V V T^ = 0. (1) 

However, as stressed by Rastall [l[ and others (see e.g. the local conservation of the 
energy-momentum tensor expressed by this relation has not been specifically tested by ob- 
servation. This is why Rastall 0] generalized the Einstein field equations by relaxing this 
assumption. Later, Al-Rawaf and Taha 7] found an alternative, equivalent form of Rastall's 
equations introducing a free parameter a and their modified equations used throughout this 
paper reduce to the Einstein ones when setting a equal to 1. Rastall's equations reduce 
to those of general relativity for empty space independently of the value of the parameter 
a. Thus, all the crucial tests of general relativity (the perihelion advance of planets in our 
solar system, the deflection of light, the gravitational red shift and the delay of radar echos) 
remain valid also in Rastall's theory. 

Rastall's modified equations becomes important mainly when treating general relativistic 
problems including sources. Since the source terms in Rastall's modified gravity equations 
include the components of the common energy-momentum tensor T^ v as well as additional 
terms of the form (1 — a) g /JU/ T/2, the source terms of Rastall's equations depend on the 
metric determined by the mass and momentum distributions of the external space. This is 
consistent with an important requirement following from Mach's principle 

In the Newtonian limit for a^l, Rastall's equations with source terms lead to a modified 
Poisson equation whose source depends on the gravitational potential, which also puts a 
strong constraint on a. It should be approximately equal to 1 in order to provide the Poisson 
equation in common weak gravitational fields. However, in large mass concentration and in 
strong gravitational field, i.e. in the stellar dynamics, even this small deviation of a from 1 
might play a significant role. 

In what follows we try to show that solutions of Rastall's equations for some model grav- 
itational systems satisfy at least one requirement following from Mach's principle, namely 
that the source terms in field equations depend, either on metric elements in case of general 
relativistic problems or on the gravitational potential in the classical Newtonian case. In this 
sense they seem to be more 'Machian' than their standard general relativistic counterparts. 

The article is organized as follows. In Section |H] we briefly summarize Rastall's theory 
rewriting his equations in the form in which the source terms are written as sums of ordinary 
energy-momentum tensor T^ v and an additional terms of the form (1 — a) g^T/2. T = T M , 
a being a free parameter. In Section IIHI we solve Rastall's equations for the isotropic 
perfect fluid calculating the components of metric tensor for arbitrary a. In Section IIVI 
we discuss the Machian property of Rastall's equations pointing out to the fact that in the 
classical "Newtonian" limit the source term (i.e. mass density) necessarily depends on the 
gravitational potential provided a^l. 

Throughout the text geometrized units in which c = 1, G = 1 are used. The metric 
signature — h ++, notation and sign conventions follow 



II. RASTALL'S GRAVITY EQUATIONS 



The standard Einstein field equations read as (see, e.g. |9|) 



Giiv = R^u — - g^uR = SnT^, (2) 
where R = -R M M ; an alternative form of Eq. (J2J) is 

R^u = 8vr (t^ - i g^T^j . (3) 

As it was demonstrated by Al-Rawaf and Taha jjj, the Rastall's modification of Einstein 
theory can be obtained by substituting a more general form of the Einstein tensor 

= aR^ v + fig^R 

for the left-hand side of the Eq. (j2J), where a and (3 are constants. In standard general 
relativity the corresponding values a = 1 and (3 = 1/2 are determined by the requirements 
that [10] 

(i) Eq. (J2J) becomes the classical Poisson equation for gravitational potential in the New- 
tonian limit and that 

(ii) energy-momentum is locally conserved, i.e. Eq. (pQ) holds. 

However, if we relax the latter condition, then the constants a and (3 are linked only 
through the relation 0] 

a (a — 2) ,3 

^W^Y a *°-2- (4) 



Using Eq. one obtains the generalized field equations first published by Rastall [1| 

Rp,v - 2 ignvR = kT^, (5) 

where 7 and k are the following functions of the parameter a 

2- a , 8vr 

7 = , k = — . 

3 — 2a a 

Evidently, if a = 1 (i.e. also 7 = 1), then Eq. © reduces to Eq. (J2J), so the Rastall's theory 
includes standard general relativity as its special case. 

The conservation condition (0) is not generally satisfied in Rastall's theory and according 
to Eq. it is replaced by [j| 

V,T" = , ) ~ a - d u R = -^—^ d u T. 

2k (3- 2a) M 2 

To make use of already known solutions of general relativity it is convenient to rewrite 
Eq. © into the conventional form employing the standard Einstein tensor 

G ^y RfJ,y ~ g^uR kT^y T^y^ (6) 
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TABLE I: The functions 7£ and "P for various values of the parameter a; let us remind that the 
values a = 0,3/2 are excluded by Eq. Q 



where 

^ = T^v — — g^T. (7) 

Therefore, it is formally possible to solve Eq. (jUJ) like the standard Einstein equations (0) 
with the modified energy-momentum tensor 7^„ which itself need not correspond to any 
physically realistic situation. For example, for the perfect fluid, the model important both 
in astrophysics and cosmology, the energy momentum tensor reads as 

Tfj,u = (q + P) u^Uu + pg^, (8) 

with the density g = g(t), pressure p = p(t) and four- velocity u fulfilling the condition 
u^u^ = — 1. Using Eq. and the fact that T = 3p — g, we can introduce the tensor 

% v = (TZ + V) UfjU u + V 9fXV , (9) 

where 

ft = i(3-a)0-|(l-a)p, (10) 

p=-I(l-a)e+i(5-3a)p. (11) 

The linear dependence of the functions 71 and V on g, p for some values of a is given 
in Table HJ it will be demonstrated below that the value a = 3 is not acceptable for our 
observed physical universe. Moreover, the Bianchi identities and Eq. (JHJ) ensure the vanishing 
divergence of introduced tensor T^, i.e. 

V V T» U = 0. (12) 

Therefore perfect fluid solutions of Rastall's equations (0) or © are formally equivalent to 
the solutions of Einstein equations (j2j) for the fluid with "density" 7Z and "pressure" V given 
by (fTUJ) and (fTTjl respectively. We exploit this fact when studying particular cases. 

In the following sections we assume that a is a universal parameter with the same value 
for all physical systems. It is always important to check that obtained results reduce to 
known solutions of standard general relativity for a = 1. 



III. STATIC SPHERICAL SYMMETRIC CASE 

First, we find the static isotropic metric fulfilling the Rastall's equations © assuming 
that the energy-momentum tensor T^ v corresponds to the perfect fluid described by Eq. ©. 
Formally, we will follow calculations analogical to those presented e.g. in j9L HlT|. 



The metric of a static, spherically symmetric spacetime can be written in the form 



ds 2 = - exp [2B (r, a)} dt 2 + exp [2 A (r, a)} dr 2 + r 2 (ch? 2 + sin 2 Mcp 2 ) 



(13) 



where we introduce two functions A (r, a) a B (r, a) and employ standard spherical coordi- 
nates of a distant observer (we implicitly suppose the asymptotic flatness of the spacetime). 
Then, the non-null components of the standard Einstein tensor read as 
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(14) 
(15) 

(16) 
(17) 



For a static fluid it also holds it* = exp(— B), ut = exp(B) and u r - u" u r 



• '' - «iV = o, therefore 

according to Eq. (JBJ) we obtain non-vanishing components of energy-momentum tensor 

„2 rri „• 2 . 



T u = gexp (2B) , T„ = pexp (2 A) , T m = pr , T w = sin i9T, 



oo 



(18) 



which can be substituted into Rastall's equations (JBJ). 

Is is convenient to solve the (t, t) component of these equations through introducing a 
different unknown function M. (r) defined as 



2M 



M(r) = - [1 - exp (-2A)] 
or compared with the metric (jl3)l 

g rr = exp (2A) = ^1 - 

Then the (t, t) Rastall equation implies 

dM k 2 4vr 2 

= r Tl = — r K; 
dr 2 a 

substituting for TZ and integrating with the condition M. (0) = 0, one gets 
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(22) 
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Evidently, setting a = 1 we come to the familiar result of general relativity 

R 



M{R)=m (R) = J Ang r 2 dr, 



where m(R) is often called as mass function and outside the fluid, i.e. for a distant observer, 
it determines the mass of the studied spherical symmetric body with a surface at r = R. 
Using this result we can transform the (r, r) component of Eq. © into the form 

dB _ M + 4vrPr 3 /a 
dr ~~ r (r - 2M) 

As we can see, the two remaining Rastall's equations for i? and (<p, if) components are 
completely equivalent. Instead of solving them explicitly, commonly another equation is 
employed instead; substituting from Eq. (jl8)l into formal "conservation law" (fT2*|) one gets 



dr dr ' dr dr 



1 . . 1 . 
-- (1 -a)g + - (5 - 3a) p 



(24) 



Again, setting a = 1 one arrives at the general relativistic counterpart of the above equation 

, dB dp 
dr dr 

Moreover, combining Eq. (}2*3*j) and (12411 we obtain the Rastall's counterpart of the 
Oppenheimer-Volkov equation (see e.g. ll2T]) 

dP = _ P) ^w. (25) 

dr r (r - 2A4) 

Here, p is always understood to be related to g by the equation of state, which also determines 
the relation between V and 1Z. In full analogy to general relativity equations (jUj) , (J2HJ) 
and (|25jl may be collectively referred to as the equations of structure for spherical stars in 
Rastall's generalization of the Einstein gravitational theory. 

Let us look for the counterpart of the famous Schwarzschild solution. Let p, g ^ only 
for r < R, where R corresponds to a surface of a spherical star. Then Eq. gives 

R R R 

MiR) = — [ TZr 2 dr = [ 47rgr 2 dr + - / Atc (g - p) r 2 dr, 

a J J 2 a J 

oo o 

where the first part corresponds to the standard Schwarzschild mass M sc h w and the second 
represents "Rastall's" correction. For a distant observer at r > R the value M. = M. (R) is 
constant and the metric coefficients read as 

g rr = exp [2 A (r, a)) = - - ^j^r , 9tt = - exp [2B (r, a)] = — M — 
where 

R 



3 X oi C 
M = M schw + - 47c{g- p) r 2 dr. 

2 a J 



o 



Naturally, putting a = 1 we get the standard Schwarzschild solution. Considering only dust 
(p = 0) solution, which is theoretically interesting, but physically highly non-realistic case, 
one come to the result 
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3 1 — a f . , 3 — a 

M schw . 
2 a J 2a 
o 



M = M schw + / 4vr^r 2 dr = — — M £ 



This puts a strong upper bound for the value of a - evidently it must hold that a < 3 to 
ensure M. > 0. 

The inconvenience of the value a = 3 can by illustrated by the following example. Let 
us calculate the metric elements A(r, a = 3) and B(r, a = 3) for an ideal perfect fluid with 
zero pressure p = (i.e. for a dust) as a source. In this case Rastall's equations (jOJ) turn 
out to be 

GW= y (T^-g^T) (26) 

with T = 3p — g = g and functions 7£, P listed in the last row of Table HI particularly 
1Z = 3p = and V = 2p + g = g. Solving Eq. (|2~T|) we come to the condition dAi/dr = 0, 
which is identically valid for all r, thus taking into account asymptotic flatness we get M. = 0, 
A(r,a = 3) = ln(l) /2 = 0, g rr = exp (2 A) = 1, which means the space-like part of the 
metric is flat. 

On the other hand, Eq. (J25|) leads to 

dB _ An 
dr 3 ^ 

Setting the conditions of asymptotic flatness B(0) = B(oo) = again, one obtains 

oo r oo r 

f dB [An C Air f An 

/ — dr = B(oo) - B(0) = = I —gr'dr' + I — gr' dr' = / — gr' dr' - B(r) 

r 

and 

r 



B{r) = / — gr dr', g tt = - exp (2B) = - exp 
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Thus the fluid affects only the (t, t)-part of the metric, while the space remains Euclidean, 
which is physically acceptable only as an approximation of weak gravitational fields. 



IV. MACHIAN PROPERTY OF RASTALL'S FIELD EQUATIONS 

An important consequence of Mach's principle is that inertial properties of a mass body 
are determined by the distribution of mass-energy in the neighbouring space. Here, the 
problem arises how to adequately characterize the physical state of this space. It seems that 
the simplest way to do it is either to take the gravitational potential in classical gravity, or 
the metric in general relativity. That is, to suppose that the source of gravitation, either 
mass-density or the elements of energy-momentum tensor, depend either on the gravitational 
potential or on the corresponding metric, respectively. 

Let us recall Rastall's equations in the form © 

Guv = Rfiu - - g^uR = — \Tuu + n a g^T) = — [T^ u + a (a) g^T] , 
I a \ A J a, 

where a (a) = (1 — a) /2. We see that the source contains elements of metric tensor. This 
is in accord with the conclusion of Mach's principle which says that the inert mass 

of a mass body depends on the magnitude and mutual distance of the neighbouring mass 



distribution. Thus, from the point of view of Mach's principle, Rastall's field equations seem 
to be more "Machian" than those of the Einstein general relativity. 

In the Newtonian gravity limit the metric in the Cartesian coordinates can be put down 
form @, EH 

ds 2 = - (1 + 2(f)) dt 2 + (1 - 2(f)) (dx 2 + dy 2 + dz 2 ) , 

where is the classical gravitational potential (typically, = —M/r outside a spherical star 
of the mass M). Taking into account the approximation of the energy momentum tensor 



Ttt = P, T 



we derive from Eq. (JHJ) 



a \ 2 J a 



1 + i— - (1 + 2(f)) 



Apparently, setting a = 1 we get the standard Poisson equation A0 = 4ng, otherwise 

A0(f) = 4irg(r) [A + 50(f)] , (27) 

where 

3(1 — a) . „ 1 — a , , . 

A = l + — - and B = 0(f). 

2a a 

Eq. ()27j) shows that in the classical limit the source term on the right-hand side of Rastall's 
equations depends on the potential of gravitational field in a specific way, which puts a strong 
constraint on a. It should be approximately equal to 1 in order to obtain the familiar Poisson 
equation for weak gravitational fields. However, for large mass concentrations and strong 
gravitational fields Eq. (|2*7jl might play a significant role in astrophysics and cosmology 0] , 
even for small deviation of a from 1. The fact that the source term of Eq. (J27j) depends 
on the potential 0(f) might be understood as an effect consistent with Mach's principle in 
non-relativity physics because the mass density g depends on the all neighbouring mass- 
energy distribution through the scalar potential 0(f). Moreover, Eq. (J27|) as well as Eq. (jHJ) 
shows that for a ^ 1 the source terms depend, contrary to general relativity, on spacetime 
coordinates, i.e. in Rastall's theory the concept of variable mass appears. 

The concept of variable rest mass in the theory of gravity is not new. It appears in 
Dicke's reformulation of Brans-Dicke theory in which the metric obeys Einstein's equations, 
but in which rest masses vary in particular way 0| or in Bekenstein's theory of the universal 
scalar field to mention only two examples. The introduction of variable rest masses in 
the theory of gravitation is part of the general program for the possible incorporating some 
features of Mach's principle into general relativity which may be reached mainly by the 
following assumptions: (i) that particles interact with some external fields as they move on 
trajectories or by Mannheim's 0] assumption of an external scalar field, which gives 
a position dependent coupling to motion; (ii) that particles move on standard geodesies 
but the Christoffel symbols are modified because the equation of motion for the metric 
is different from the standard one. This can be reached either through a fundamental 
replacement of the Einstein curvature tensor by another equally covariant second-rank 
tensor, or through introducing a more complicated energy-momentum tensor than it is 
conventionally considered, keeping the common Einstein tensor at the same time. 



Concluding remark 



Rastall modified Einstein' equations by replacing the second term in the Einstein ten- 
sor. The advantage of this substitution is that Rastall's equations reduce to Einstein's one 
for empty space-time. Differences between Einstein's and Rastall's equations arise when 
one investigates gravitating systems with non-zero sources. We have showed that consid- 
erable differences between both theories exists also regarding the incorporation of Mach's 
principle in the theory of gravitation. Hence, Rastall's equations represent a set of gravity 
equations which includes also the Einstein ones. In this sense, they can be considered as a 
generalization of general relativity. 
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